Interaction induced ferro-electricity in the rotational states of polar molecules 
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We show that a ferro-electric quantum phase transition can be driven by the dipolar interaction of 
polar molecules in the presence a micro-wave field. The obtained ferro-electricity crucially depends 
on the harmonic confinement potential, and the resulting dipole moment persists even when the 
external field is turned off adiabatically. The transition is shown to be second order for fermions 
and for bosons of a smaller permanent dipole moment, but is first order for bosons of a larger 
moment. Our results suggest the possibility of manipulating the microscopic rotational state of 
polar molecules by tuning the trap's aspect ratio (and other mesoscopic parameters), even though 
the later's energy scale is smaller than the former's by six orders of magnitude. 
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The successful realization of high phase-space-density 
polar molecules in the JILA group 1] opens a new di- 
rection of strongly correlated quantum gases: ultracold 
polar molecules. Many experimental [2| and theoreti- 
cal works have been carried out in recent years. 
Most many-body physics proposed so far are based on 
the situation when the inter-molecule effective interac- 
tion is predetermined (by a separated calculation) Q and 
cannot be affected by the many-body physics. On the 
other hand, the reverse operation, i.e. manipulating the 
molecular rotational state by changing a mesoscopic pa- 
rameter (say trapping frequency or particle number etc.), 
is still unfeasible, because the rotational energy is about 
six orders of magnitude larger than the interaction energy 
in a dilute gas. This makes systems of polar molecules 
very different from spinor atoms Q, although they both 
have rich internal states and long-ranged dipolar inter- 
action. For the same reason, a spontaneous long-ranged 
order of molecular rotational state seems impossible ei- 
ther while the related many-body phenomena, like 
(anti-)ferro-magnetism and multi-ferroics [lOj . have been 
important subjects of solid state physics for decades. 

In this paper, we demonstrate that the rotational state 
of polar molecules can be manipulated by mesoscopic pa- 
rameters through an external AC field Q: When the 
trapping potential is elongated and the field frequency is 
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FIG. 1: (Color online) (a) The lowest two rotational eigen- 
states of a polar molecule. Here we consider a linearly po- 
larized AC field to couple state (\g) — |0, 0}) and state 
(je) = 1 1 , 0} ) with detuning A. (b) Schematic figure for the 
ferro-electricity in the presence of a microwave field. The blue 
up-down arrows indicate the oscillating dipole moments. 



tuned close to the rotational energy (see Fig. [Ha)), the 
internal rotational state of molecules can all be "locked" 
to each other by interaction, leading to a spontaneous 
mesoscopically ordered state, ferro-electricity. Different 
from the standard dressed state in a single molecule pic- 
ture Q, the interaction-induced macroscopic dipole mo- 
ment here sustains even when the external field is turned 
off adiabatically. We calculate the quantum phase dia- 
grams for both bosonic and fermionic polar molecules in 
an elongated ellipsoid trap. We then discuss the nature 
of the phase transition and how it can be tuned via meso- 
scopic parameters, like trap's aspect ratio or the detuning 
of an external field. 

The system Hamiltonian of N polar molecules in- 
cludes both the rotational and the orbital degrees of 
freedom: H to t = Yli Hrot,i + H orh . Here H Iot ,i = 
BJ 2 — di ■ E(t) describes the rotational state of the 
ith molecule with JL and B being the angular mo- 
mentum operator and the rotational constant respec- 
tively. E(t) is the external field, and dj is the elec- 
tric dipole moment. The orbital part of the Hamil- 
tonian is given by -ff rb = Y^iL 
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position and momentum operators, and m is the molecule 
mass. N is the total number of particles and Vt(r) = 
^m(ujp(x 2 + y 2 ) + cu 2 z 2 ) is the external trapping poten- 
tial. The second term is the dipolar interaction with 
e,-j = ry/|r^| being the unit vector of the relative posi- 
tion, Yij =Vi— Tj . Following the spirit of the mean-field 
approximation, without losing generality, we can approx- 
imate the ground state wavefunction by a factorizable 
form: *(1,2,---,A) = $ orb (0; n • • • r N ) ® nili l*i>i> 
where $ O rb(0; ri • • • rjv) is the orbital wavefunction with 
a variational parameter, and \0i)i describes the rota- 
tional state of the ith particle as more well-defined below. 

In this paper, we are considering the situation 
when a linearly polarized AC field is applied, E(t) = 
Eac cos(ujt)z with an angular frequency w, to couple ro- 
tational states \g) = |0, 0) and |e) = 1 1 , 0) (see Fig. QJa)). 
The single molecule state can be solved by the Floquet 
theory [12j | within the rotating wave approximation, and 
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the following time-independent Hamiltonian (in the trun- 
cated rotating basis, {\g),e~ tult \e}}) can be obtained: 
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A = 
Jl 2 , and 



where SI = d^Ep^fh is the Rabi frequency. 
uj — 2B/h is the detuning, Aq = -v/A 2 
tana = ft/ A. Here do = |(g|di|e)| ^ s the perma- 
nent dipole moment. According to the Floquet the- 
orem, a general rotational state can be expressed to 
be \0t)t = cos (U) \v+) e- iE + t ' h + sin (|) \ V -) e 



-iE-t/h 



where E± — (A ± y/ A 2 + ft 2 ) are the eigenener- 
gies, and |u + ) = — sin(a/2)|p) + cos(a/2)e _i "*|e) and 
|t!_) = cos(a/2)|g) + sin(a/2)e~ l " t |e) are the correspond- 
ing eingenstates of H rot ^. After a straightforward calcu- 
lation, we can obtain the rotational energy per parti- 
cle, E mt {9)/N = -M _ OAs. J2? =1 cos6i, and the av- 
eraged dipole moment (in the limit of small Aq/lv): 
(di(t)) = (i9i\di\9i)i — dosm(9i — a) cos(ut)z, which is 
oscillating in time. If hA^ is much larger than any other 
orbital energy scale, the rotational energy is then always 
minimized at a single dressed-state, \v+) (i.e. 9i = 0) 
with a finite rotating dipole moment. However, this state 
alone is not a true ferro-electric state @ , no matter what 
kind of orbital wavefunction it has, because the dipole 
moment disappears in the limit of zero external field (i.e. 
a — > 0). However, the situation can be totally different 
when A is tuned to be comparable to the interactional 
energy, as shown in more details below. 

Before further investigation, we note that the 9i, in 
principle, can depend on the confinement geometry: 
when the molecule cloud is elongated along the z axis (i.e. 
the polarization direction), 9i will tend to be uniform in 
space in order to gain the interaction energy. However, if 
the cloud is oblate in the x — y plane, 0i can become non- 
uniform in space, i.e. for ming a domain wall structure or 
a rotation-orbital texture [16j . Since here we concentrate 
on the ferro-electric order rather than the complete phase 
diagram, we will always consider a highly elongated trap 
and assume a uniform polarization for simplicity. A full 
theory including non-uniform polarization will be pre- 
sented in the future. 

Within such approximation, the interaction strength 
can be characterized by the following dipolar length scale 
in the dilute limit jj], 

_ m((d(t)) ■ (d(t))) t 4 0) 
a d (9) = — = — sin (9 -a), (2) 

where a d = ind^/H 2 , and (• • • ) t = ~ /J" dt ■ ■ ■ is the time 
average over many periods, r 3> 2it/uj. To study its effect 
on the orbital many-body wavefunction, we first consider 
bosonic molecules at zero temperature, where the orbital 
condensate wavefunction can be well- approximated by a 
Gaussian type form: *(r) = n 1 ' 2 -( P 2 +z 2 /p 2 )/R 2 p 

with Rp being the horizontal radius and (3 being the as- 
pect ratio. The orbital condensate energy, E^ c (6), then 
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FIG. 2: Results of bosonic polar molecules: (a) Total groud 
state energy as a function of 9. Curves from top to bottom are 
respectively for hA/xt, = 6, 5, 4, and 3. (b) The calculated 
effective dipolar length, a<i(8), as a function of the detuning, 
A. Curves are for /a 3 = 5.5, 4, and 2, respectively from 
top to bottom, (c) Quantum phase diagram as a function of 
af and A. The para- to ferro-electricity transition is second 



order for /a s < 5, while it is first order for a^'/a s > 5. 
When /a s > 6, the condensate becomes unstable toward 
collapse, (d) Critical detuning for the second order transition 
line as a function of trapping aspect ratio, n. Curves from 
top to bottom are for /a s = 5, 4, and 3. Note that we 
have taken fl — ¥ for all results shown here. 
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where a p = yjh/mujp is the oscillator length, and a s is 
the s-wave scattering length, k = lo p /uj z . 

Now we minimize the total energy, E^ c (9) = 
E TO t(6) + E^ c (9), with respect to variational param- 
eters, R p , (3, and 9. For an elongated trap (i.e. /3 3> I), 

we can use A2(P) — > and therefore the interaction 
energy term of Eq. (J3j) becomes the same as the re- 
sult of a non-dipolar condensate with an effective length, 
a e g(9) = a s — ad(9)/3. Within the Thomas-Fermi limit 
(i.e. AT 3> 1 and hence the kinetic energy is negligible), 
the minimum energy at a given 9 can be done analytically 
[Ilf , We obtain j3 = k = uj p /uj z , and 
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where Q = (w 2 ,^) 1 / 3 and a = \J h/mCo. It is easy to sec 
that the system becomes unstable toward collapse when 
a cff($) < 0, giving the same condition derived in mag- 
netic atoms [14| . Taking zero field limit and expanding 
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E^ c (9) to the leading order of small 9, we obtain 
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after shifting the total energy to zero at 9 = 0. Here 

(0) 

we define Xb = j^rXb to be the "electric susceptibil- 

ity" and A = ± - Here X ° = ^(O)^" 1 . It is 

easy to see that a second order quantum phase transition 
from a para-electricity toward a ferro-electricity can take 
place when |A| < A c = ixb/h and when \\b — ^§f- > 

(i.e. a d 0) < 5o a ). In Fig. I^a), we show E™ c (9) as a 
function of the variational parameter, 9, near the phase 
transition of of /a s — 2. However, when > 5a s , 
the coefficient of the 9 A term becomes negative, lead- 
ing to a first order phase transition instead, see Fig. 
[2]Jb). In Fig. HJc), we show the calculated quantum 
phase diagram for such elongated condensate, as a func- 



tion of a d °^ and A. A critical point is at a d "' ja s = 5 

and hA/x® ~ 1-33. When a d /a s > 6, the conden- 
sate becomes collapse due to the strong attractive dipo- 
lar interaction. Note that since E^ c (0) = N\° de- 
pends on the trap's aspect ratio, k = uj p /uj z , we can 
also derive the critical detuning as a function of k to be: 

hA c = ^jXb = ^t^b( K = 1)^ 2/5 , and therefore the 
ferro-electricity can be also tuned easily by changing the 
trap's aspect ratio as shown in Fig. ^d). Note that 
the decay of A c for larger k is simply due to the reduced 
particle density. Taking 39 K 87 Rb as an example, we have 

do = 0.5 Debye and hence a d °^ = 0.47 /im. Using a s = 0.2 
fxm, we can easily obtain A c = 2ir x 250 Hz for N = 10 4 , 
Lo p = 2tt x 100 Hz, and k = 100. 

Now we turn to systems of fermionic molecules. To 
the best of our knowledge, so far only a variational study 
[l^ | and Fermi-liquid theory in a uniform system [Tsj] 
are studied in the literature, while a full self-consistent 
Hartree-Fock theory with the local density approxima- 
tion (SCHF+LDA) have not been investigated in the lit- 
erature. It is therefore worthy to show our results on this 
aspect first. Within such HF+LDA approximation, the 
orbital energy can be calculated to be 
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where S(k,r) = J r , k , \ V d {v - r') - V d (k - k')5(r - r' 
/k'(r')is the static local self-energy, and /k(r) = 
[exp ((e° + V t {r) + S(k, r) - n)/k B T) + l] ^ 
distribution at temperature T. J r . = 
is the integration over the whole phase space. 
Vrf(r) = |(d)| 2 (r 2 — 3z 2 )|r|~ 5 is dipolar interaction 
and Vd(k) is its Fourier transform in momentum space. 
e k = ^- 2 k 2 /2rn is the kinetic energy, and /x is the chemical 
potential. Eq. © can be also applied to systems of 
finite temperature, but in this paper we just consider 
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FIG. 3: (a) Zero temperature density profile (after column 
integration) for fermionic polar molecules, 40 K 87 Rb, in three 
different trap's aspect ratio, ks. Here we use JV = 3 x 10°, 
|(d) | = 0.5 Debye, and a fixed lo p — 2ty x 100 Hz. (b) and 
(c) are renormalized cloud's aspect ratio in the spatial and 
momentum density distribution as a function k. |(d)| = 0.1, 
0.3 and 0.4 Debye from bottom to top. 



zero temperature situation. In Fig. pj^a), we show the 
calculated density distribution for different aspect ratios 
of the trap, and in (b) and (c) we show the obtained 
cloud's aspect ratio (normalized by the non-interacting 

,!/2 

and 



aspect ratio). Here R z = (± J r k /k(i")^ 

K z = I -k J , /k( r )^ 2 ) are the averaged radii in the 

spatial and momentum space (similar definitions for 
R x and K x ) . We find that our results are consistent 
with Ref. [l9| for spatial aspect ratio, but disagree with 
them in the momentum aspect ratio. It is probably 
because their variational approach did not capture the 
full contribution of Hartree energy. 

From Eq. ([5]), one can see that the susceptibility has 
two contributions: one is from the self-energy in Eq. 

and the other is from the single particle energy via 
the deformation of Fermi surface. In low temperature 
limit, the latter contribution can be neglected, because 
it involves only an integral on the Fermi surface of the 
phase space. Besides, most effects of the self-energy in 
/k(r) can be absorbed by the shift of chemical potential 
due to the conservation of total number of particles. As 
a result, we can just keep the interaction effect to the 
leading order of aa(9), which is to replace m\(d)\ 2 /h in 
Vd(r) of Eq. ©. The first order transition does not 
exist here for the similar reason, since the high order 
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FIG. 4: Results of fermionic polar molecules: (a) Quantum 
phase diagram in terms of trap's aspect ratio, k, and the 
detuning A, in the zero field limit, (b) The calculated dipolar 
length, a,d(0), as a function of A. Results for different Rabi 
frequencies are shown together. 



behavior is enhanced more strongly in a bosonic gas. 
After re-scaling the length scale to rp = J 2e F °/mui 2 

and the momentum scale to kp = ^j2me F °/h 2 with 
(67V) 1 / 3 hu> [H[, we obtain (taking O ->■ 0), 
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where 

is a dimensionless function of k with Yi m (e r ^) being 
the spherical harmonics, efc is the unit along k — k'. 
J k stands for the integral inside an elliptic sphere, 

x 2 +y 2 + z 2 /k 2 + k 2 = 1. Note that C(n) plays a similar 
role as A2(j3) in Eq. Therefore, expanding the 

total energy to the quadratic order of 9, we can obtain 
the "susceptibility" per particle to be: Xf = XfC{ K )> 
where x°f = 48(3Y) 1/2 d 2 ,/a 3 . The critical detuning 
for ferro-electricity is given by A c = Axf/h as shown 



in Eq. ©. For 40 K 87 Rb molecules with N = 10 4 , 
Lo p = 2tt x 100Hz, d = 0.5 Debye and n — 100, we 
have x°f = 44.4 kHz, and C(k) = 3x 10~ 4 and hence 
A c = Axf/fr ~ 53 Hz. This value is smaller than the 
corresponding value of bosonic molecules, because the 
condensate enhances the dipolar interaction effects. 
In Figs. [4ja) we show the calculated quantum phase 
diagram, which is qualitatively similar to the bosonic 
case (see Fig. [Did)). In Fig. Q|b), we show the dipolar 
strength as a function of detuning. Results for different 
Rabi frequencies are also shown together. 

Finally, we note that the ferro-electric order predicted 
here can be easily observed from the electric field dis- 
tribution when the AC field is turned off adiabatically 
(a —¥ 0). The resulting macroscopic polarization can also 
be manipulated by changing the confinement's aspect ra- 
tio, s-wave scattering length (bosonic case), or detun- 
ing frequency. Such macroscopically coherent rotational 
state can be neither reproduced in a single molecule pic- 
ture nor in a uniform system. Finally, the "spin" -wave 
fluctuation as well as the texture in the rotational state 
shall also be an interesting subject for future study. 

In conclusion, we have shown that an interaction 
driven quantum phase transition from a para-electric to 
a ferro-electric quantum gas of polar molecules can be re- 
alized in the presence of a micro- wave field. We study the 
nature of such quantum phase transition, and then pro- 
pose how to manipulate the molecular rotational states 
by the tuning of many-body parameters. 
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